The low-frequency response in the surface superconducting state of ZrBi 2 single 

crystal 
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The large nonlinear response of a single crystal ZrBi2 to an ac field (frequency 40 - 2500 Hz) for 
-ffo > Hc2 has been observed. Direct measurements of the ac wave form and the exact numerical 
solution of the Ginzburg-Landau equations, as well as phenomenological relaxation equation, permit 
the study of the surface superconducting states dynamics. It is shown, that the low frequency 
response is defined by transitions between the metastable superconducting states under the action 
of an ac field. The relaxation rate which determines such transitions dynamics, is found. 

PACS numbers: 74.25.Nf; 74.60.Ec 



I. INTRODUCTION 



Recently high-quality superconducting ZrBi2 single 
crystals with transition temperatures T c — 6.06 K have 
been grown. The investigation of their physical proper- 
ties, including electron transport, tunnel characteristics, 
and critical fields, have shown that the Ginzburg-Landau 
parameter re is only slightly larger than the boundary be- 
tween the type-I - type-II superconductor value 0, 0- 
In this paper we concern the low-frequency response 
of a ZrB 12 crystal when the dc external magnetic field 
Ho > H C 2 is parallel to the sample surface. In spite 
of the fact that the sample is in the surface supercon- 
ducting state (SSS) 3!], no static magnetic moment is 
observed, while the ac response in this regime is large 
and nonlinear even for an ac amplitude ho « Ho- In- 
deed, at equilibrium the total surface current equals zero 
in SSS, the internal dc magnetic field in the bulk equals 
Ho, and the magnetic moment of the specimen is small. 
On the other hand, the ac magnetic field drives the sam- 
ple into a metastable SSS where the total surface cur- 
rent reaches a finite value. The internal magnetic field 
deviates from the external one, and as a results the ac 
response becomes large. The low-frequency response of 
superconductors in SSS was the focus of intensive exper- 
imental investigations jj, la, |6J since the first prediction 
of the existence of SSS in The observed wave 
form of the ac response, corresponding to the flux pass- 
ing through the specimen, explicitly invoked a model sim- 
ilar to the Bean model Q- Recently SSS attracted re- 
newed interest from various directions as described in 
Ref. [1,111 El El El El- Paramagnetic effect in a super- 
conducting disk 0, stochastic resonance 0, , the per- 
colation transition in the field Ho = 0.81-ff C 3 El have 
been observed. It was proposed to use low-frequency 
response for testing the quality of superconducting res- 
onators in accelerators |ll| . Surface states were observed 
also in single crystals of MgB2 Our experimental 

results presented here show that in ZrBi2 single crystal, 



the Bean critical model of surface sheath does not give 
an adequate description of the observed wave form which 
corresponds to the flux passing through the sample. In 
the framework of Ginzburg-Landau, theory we calculated 
the surface current in metastable SSS's which exists un- 
der an ac magnetic field. Observing the wave forms, we 
studied the metastable SSS dynamics and determined the 
relaxation rate under an ac field. We found that the re- 
laxation time for transition to the equilibrium state is not 
constant and depends on the surplus of the free energy. 
This relaxation time is decreased with the dc magnetic 
field, and depends on the driving field frequency. 



II. EXPERIMENTAL 

The measurements were carried out at T = 5 K on 
ZrBi2 single crystal. The sample was grown in the In- 
stitute for Problems of Materials Science, Ukraine. Its 
dimensions are 10.3 x 3.2 x 1.2 mm 3 and it was cut by 
an electric spark from a large crystal of 6 mm diame- 
ter and 40 mm length. The surface of the sample was 
polished mechanically and then chemical etched in boil- 
ing HNO3/H2O (1:1) for 10 minutes was used. X-ray 
pictures showed that a sample was single-phase material 
with the UB12 structure (space group Fm3m, a = 7.407 
A [lj|. The tunnel characteristics of this sample were 
described earlier ■ The dc-magnetic moment was mea- 
sured using a SQUID magnetometer. A block diagram 
of the ac linear and nonlinear setup is shown in Fig. ^ 

The ac magnetic field h(t) — ho sin(wi) was supplied 
by the magnetometer copper solenoid. The ac response 
was measured by an inductive pick-up coil method [l5| . 
The sample was put into one coil of a balanced pair of 
pick-up coils and the induced voltage V(t) oc dM(t)/dt 
was measured with an oscilloscope. Here M is the mag- 
netic moment of the sample. The lock-in amplifier was 
used in order to measure simultaneously in-phase and 
out-of-phase signals of the first and third harmonics of 
the driving frequency. An oscilloscope measured the wave 
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FIG. 1: (Color online) Block diagram of the experimental 
setup. LFG - low frequency generator, LCK - lock-in amplifier 
and OSC - oscillograph. 



form of the signal in one channel. The second channel 
of the oscilloscope measured the time derivative of the 
excitation field. 



0.01 




III. EXPERIMENTAL RESULTS 



Fig. shows the real part of the ac susceptibility at 
the fundamental frequency, x > an d the zero-field cooled 
(ZFC) dc susceptibility, Xdc — M/Hq as a function of dc 
field Hq. The inset to Fig. presents the ZFC mag- 
netization curve at 5 K. Field dependencies of the ac 
susceptibility imaginary part at fundamental frequency, 
X , and amplitude of the third harmonic, A 3uj are shown 
on Fig. [2)3 and Fig. |2t respectively. Amplitude depen- 
dence of the third harmonic, A 3LJ (ho) for Hq = 180 Oe 
is presented on inset to Fig. It is clear that this de- 
pendence is far to be cubic as the perturbation theory 
predicts. Experiment shows that amplitude dependence 
of A 3LU (h ) is not cubic at any dc field Hq. 

It is clear that the observed large signal of the A 3uJ 
and maximum of the x located in a magnetic field 
H C 2 < Hq < H c3 , e.g. in a surface superconducting 
state, although the zero dc signal indicates that bulk of 
the sample is in the normal state. The absorption in the 
SSS exceeds losses in the mixed and normal states. 

Fig. |2| shows the time derivative of the magnetic mo- 
ment of a sample at different applied magnetic fields at 
T = 5 K. Note (i) that only in the SSS, the signal does 
not have the sine-form, (ii) The amplitude dependence 
of the third harmonic, A^iho), presented in the inset to 
Fig. |2t does not exhibit any cubic dependence. 

The experimental data presented in Figs. |21 and [3] are 
complex and the theoretical model which explains these 
observations is given in the next section. 
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FIG. 2: (Color online) (a) \ and Xdc = M/Ho magnetic field 
dependencies at T = 5 K. Inset: magnetization curve after 
ZFC. 

(b) Magnetic field dependence of x" ■ 

(c) Field dependence of the amplitude of the third har- 
monic, A^. Inset: amplitude dependence of Azwiho) at 
Ho = 180 Oe. 

ac measurements were carried out at frequency lu/2tt = 
170 Hz and ho = 0.4 Oe. 
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FIG. 3: (Color online) Oscillogram of dM/dt for three differ- 
ent magnetic fields in the Meissner state (Ho — 0) , in surface 
superconducting state, (-ffo = 180 Oe < H C 3), and for a nor- 
mal state (H = 300 Oe > H c3 ) at ho = 4.75 Oe. 



IV. THEORETICAL MODEL 

Our theoretical approach is based on the numerical 
solution of the two Ginzburg-Landau equations f° r 
the order parameter and vector potential, which in the 
normalized form are as follows: 



-(zV/re- 
-curlcurl^l = 
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The order parameter, 'J', is normalized with respect 
to its value in zero magnetic field, the distances with re- 
spect to the London penetration length A, and the vector 
potential A with respect to ^/2H C X, where H c is the ther- 
modynamic critical field and k — A/£ is the Ginzburg- 
Landau parameter, and £ is the correlation length. In 
the second equation for the vector potential we neglected 
the normal current, assuming that the skin depth ex- 
ceeds both the London penetration length and sample 
thickness. It is assumed that the sample form is a slab 
with 2d thickness and that the external magnetic field is 
parallel to its surface. The chosen coordinates are: the 
x-axis normal to the slab, (thus the symmetry plane is 
x = 0), and the z-axis is directed along magnetic field. 
It is assumed that external magnetic field H > H C 2 = K. 

Assuming the surface solutions have the form of 



^(x,y) = f(x)exp(inky) 
therefore, Eqs.^ reduce to 

f (A - kff -f + f 
B = f 2 (A~kl 



ld 2 f 

K dx 2 



(2) 



(3) 



where k is constant. The boundary conditions at x = 
±d are: df(±d)/dx = 0, dA(±d)/dx = H ; at x = 
/(0) = and A(0)=0; H = H +h sm(u;t). An additional 
condition for the surface states df(0)/dx = is satisfied 
only asymptotically for d — ► oo. For the equilibrium 
state the value of k = k eq can be obtained by minimizing 



the Gibbs free energy defined as: 

F = jdV{\\^\ A - |*| 2 + |*V*//s + A* | 
+B 2 - 2BH}, 



(4) 



where B is the magnetic induction. Using Eq.[3]one then 
obtains 



F = -HA(d) - J d dx{\f\x) + A(x)[A(x) 
-k]p(x)} 



(5) 



The two coupled Eqs.|3] could be solved by numerical 
methods. The order parameter for surface solutions de- 
viates from zero only near the sample boundary, and we 
could consider comparatively small d/X < 10. The actual 
sample thickness exceeds A by 4 or 5 orders of magnitude. 
The solutions for large d could be found from the ones 
for small d by transformation k = k s + Hi(d — d s ), where 
Hi = H s (0, H, K) is the magnetic field at x = in the 
problem for d = d s = 10A. The index s corresponds to 
the solution for this small d. This choice of d s is sufficient 
for numerical calculations and provides the solutions with 
/s(0) = 0, df s (0)/dx = 0. The free energy transformed 
as 



F = F S - Hi(d - d s )(H - J [A s (x) - k s ]f s (x)dx) (6) 

To simplify the calculations we use below variable k s 
and omit index s. The properties of the equilibrium 
solutions for a semi-infinite half-space have been dis- 
cussed in 16] . The order parameter, the supercurrent, 
and the internal magnetic field were calculated. In these 
states the total surface current equal zero and free energy 
reaches a minimum value. The ac magnetic field drives 
the superconductor into a metastable state. These states 
correspond to the solutions of Eqs.[3]for k ^ k eq . The so- 
lution of Eqs.|3] shows that surface states exist in a wide 
range of k near k eq as shown in the upper panel of Fig. 
but only for \k— k eq \ << k eq free energy of these states is 
lower than the energy of the normal state. Moreover, this 
range shrinks with increasing sample thickness (Fig.[5j). 

This is due to the increase of the contribution of the 
first term in Eq. which is the order of the Gibbs en- 
ergy of the normal state Fq — —H 2 d. The total surface 
current equals zero for equilibrium k and increases with 
increasing \k — k eq \ as is shown in the lower panel of 
Fig. For a given k, the free energy versus magnetic 
field does not exhibit any minimum in the equilibrium. 
Only the difference between the Gibbs energy of the su- 
perconducting and normal states exhibits minimum near 
equilibrium field as was discussed in but for a such 
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equation that governs the dynamics is 
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FIG. 4: Upper panel. The maximal value of the order pa- 
rameter, fmax, as a function of k. Lower panel. Total surface 
current, J, as a function of the k parameter. The equilibrium 
value of k = 8.8 corresponds to zero surface current. 
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FIG. 5: (Color online) Free energy of the surface supercon- 
ducting states relative to the free energy of normal state (Fo) 
in magnetic field Ho = 1 and k = 0.75 for two values of the 
sample thickness. 



representation the reference point moved with changing 
field. 

The magnetic moment of the sample actually depends 
on the total surface current J, because the current is 
localized in a thin surface layer. This current is a function 
of the external magnetic field H and k, J — J(H, k). The 
response of the sample to the ac magnetic field depends 
on the dynamics of k. A priori, one can assume that the 



dk 
~dt 



= -u[k-k eq {H% 



(7) 



where k eq {H) is the equilibrium k in instantaneous 
value of external magnetic field and v — u{k— k eq ) is the 
relaxation rate. Function k eq (H) has to be found from 
Eqs. El and for \H — H Q \ « H is well approximated by a 
polynomial of the third order of h = H — Hq. Using the 
function J(h, k) calculated from Eqs. [3] and Eq. one 
can obtain the time evolution of the surface current in 
an ac field and compare with observed wave forms. The 
time derivative of the surface current is proportional to 
the observed signal V. The coefficient a = depends 
on the experimental apparatus parameters. Actually we 
could obtain k(t) directly from experimental data and 
test the correctness of Eq. [7| We may write 



dJ(h,k)dh dJ(h,k)dk 



dh dt 



dk dt 



aV{t) 



(8) 



This expression permits us to obtain k(t) from the ob- 
served wave form. It is a first order differential equation 
for k(t). To evaluate k(t) we have to know k at t = 
and a, since the derivatives dJ/dh and dJ/dk can be cal- 
culated from Eqs. |3J The fc(0) value can be found from 
the condition when the maximal current value during the 
period is minimal. In order to find a we calculated J(t) 
assuming that v in Eq. 0is constant. Then, we choose 
the v value in order to minimize difference between the 
calculated and experimental data. This procedure gives 
both v and a. To be sure that v is actually constant, 
one has to collect the weak ac field data. The observed 
signal during one period of the ac field and the result of 
a simulation with Eqs. and [7| are shown in Fig.0] The 
data in this figure were collected in a dc field of 130 Oe, 
and an ac field with amplitude 1.78 Oe and frequency 
uj/2tt = 733 Hz. In our calculations we took v/uj = 0.05 
and the Ginzburg-Landau parameter n = 0.75 2J. The 
good correlation between the calculated and experimen- 
tal data permits one to find the scale coefficient a which 
is used below. 



V. DISCUSSION 

As was shown above the losses are small both in 
the mixed and normal states and have a maximum at 
H > H c2 (see Fig. |2t). The H c2 = 126 Oe is deter- 
mined from the dc magnetization curve (inset to Fig-EJi). 
The oscillogram, Fig. |3 in both the Meissner and nor- 
mal states (Ho = and 300 Oe) has a sine shape, and 
for H C 2 < Hq < H c s the wave form deviates from a sine 
shape. We do not observe any clear plateau for dM/dt 
in the ac period. Such a plateau is a peculiarity of the 
Bean model when it is applied to surface currents 0, Il7| . 
Using the experimental data and the model developed 
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FIG. 6: (Color online) The time derivative dk/dt plotted as 
a function of k — k eq for at Ho = 130 Oe and lo/2-k — 733 Hz 
(a): ho = 1-78 Oe and (b): ho = 5.9 Oe. The linear fit of 
the experimental curve (a) gives v/u = 0.051. The hysteresis 
at large amplitudes (b) shows that, generally speaking, dk/dt 
depends on k and the instantaneous magnetic field not only 
through k = k eq . 



in previous section one can calculate dk / dt as a function 
of k — k eg . Fig. E| shows dk/dt plotted as a function of 
k — k eq obtained from the wave forms that were observed 
for H = 130 Oe and u/2tt = 733 Hz. The linear fit 
of dk/dt at h = 1.78 Oe yields u(Q)/u = 0.051 which 
agrees well with the v/oj — 0.05 used in Eq. \7\ when 
the scale coefficient has been found. Visible hysteresis in 
Fig. 03 indicates that at a high amplitude of excitation 
the relaxation rate in Eq. [7| v depends on k and on the 
instantaneous value of h(t) not only through k — k eq (h). 

The expression for v(k — k eq ) could be found from fit- 
ting of dk/dt by the polynomial of k — k eq . The approx- 
imation expression which has a form 



FIG. 7: (Color online) The observed and calculated (solid 
lines) oscillogram for Ho = 130 Oe, lu/2tt — 733 Hz at ho = 
1.78 Oe and h = 5.9 Oe. 
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u(x)/u = 0.051- 0.1 17x+ , , 

+ 1.323.t 2 + 0.184x 3 - 0.747x 4 ^ ' 

provides the calculated data, which with an accuracy 
of better than 10%, reproduces the experimental data for 
a dc field of 130 Oe and a frequency of 733 Hz as is shown 
in Fig. E| 

Increasing the dc field leads to the increasing of the 
relaxation rate v. We found that for Ho = 138 and 180 
Oe, the relaxation parameter for weak ac amplitudes are 
v(0)/uj = 0.145 and 4.725 respectively (see Fig. El. 

The calculated wave forms with the help of the pro- 
posed model reproduce experimental data only for a weak 
ac field as shown in Fig. ED This is due to the increase 
in the difference between the two values of dk/dt for the 
same k — k eq at larger ac amplitudes, ( see at Fig. [HJj, d). 



FIG. 8: (Color online) dk/dt as a function oik — k eq at u)/2ir = 
733 Hz. (a) - H = 138 Oe, ho = 0.59 Oe; (b) - H = 138 Oe, 
h = 5.9 Oe; (c) - Ho = 180 Oe, ho = 0.59 Oe; (d) - H = 
180 Oe, ho — 5.9 Oe. Linear fit at low amplitude of excitation 
( (a) and (c) panels) gives v/u = 0.144 for Ho = 138 Oe and 
v/u = 4.73 for H = 180 Oe. 



Fig. 1 101 shows the dk/dt(k — k eq ) dependence for differ- 
ent frequencies at Ho — 130 Oe and ho — 4.75 Oe. One 
may conclude from this figure that the relaxation rate 
v (if Eq. [7| could be applied) increases with excitation 
frequency lu. 

It is clear that the model equation where the relax- 
ation rate v depends on the one variable k — k eq , is valid 
only for small ac amplitudes. The transition from the 
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FIG. 9: Calculated (solid lines) and experimental (triangle) 
oscillograms at ui/2n = 733 Hz. (a) - Ho = 138 Oe, ho = 
0.59 Oe; (b) - Ho = 138 Oe, h = 5.9 Oe; (c) - H = 180 Oe, 
h = 0.59 Oe; (d) - H = 180 Oe, ho = 5.9 Oe. 



ac amplitude and/or the dc field results in the explicit 
dependence on both: k and the instantaneous value of 
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FIG. 11: Calculated dependence of the Gibbs energy as a 
function of the k — k eq parameter during the ac cycle for dif- 
ferent values of the dc field Ho and amplitude of excitation, 
ho at frequency lu/2-k — 733 Hz. 
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FIG. 10: (Color online) dk/dt as a function of k — k eq at 
Ho = 130 Oe and ho = 4.75 Oe for different frequencies. 



surface state with one k to another fc, requires changing 
the order parameter in the whole sample. Possibly, it 
happens through the nucleation of a new phase. This 
process is governed by excess energy. In general the re- 
laxation constant in Eq. \7\ may depend on the energy of 
the state. Approximately it could be taken into account 
by the assumption that the relaxation constant depends 
on k—k eq . We see that for small ac fields and not far away 
from H C 2, this assumption is correct. But increasing the 



the magnetic field Hq. The straightforward calculations 
of the Gibbs energy F, exhibited in Fig. El shows that 
when this energy is a single- valued function of k — k eq , 
the simulation with Eq.[7| gives acceptable results. 

When F becomes a multivalued function, the calcu- 
lated wave forms differ from the experimental data. In 
order to obtain a proper theoretical description, one has 
to take into account that the energy of SSS is not ex- 
pressed only through k — k eq . 



VI. CONCLUSION 

We experimentally investigated the dynamics of the 
surface metastable superconducting states of ZrBi2 in ac 
fields at low frequencies (40 - 2500 Hz). It was shown that 
for low ac amplitudes of excitation these dynamics are 
governed by a simple relaxation equation. The relaxation 
rate depends on the deviation from the equilibrium state. 
Decreasing of the frequency of the applied ac field results 
in increasing the relaxation time. 
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